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Abstract: We give results for the distribution and number of flux vacua of various types, 
supersymmetric and nonsuper symmetric, in lib string theory compactified on Calabi-Yau 
manifolds. We compare this with related problems such as counting attractor points. 



Contents 



0- 



Introduction 
|l.l| Distributions of vacua 

Approximate distributions of vacua 
Other types of vacuum 
Index densities 
Ensembles of flux vacua 



1.4 



1.5 



|2|. Notat ions and some useful formulas 

.1| Orientifold limit 



|3|. Distributions of supersymmetric vacua 

Computing densities 
|3.1.l| The case n = 1 
Index density 



3T| 



3.1.4 



3.1.5 



Example 1 
Example 2 
Example 3 



Conifold 
Mirror Quintic 



3.1.(: d |VF| signature distribution 



3.2 



3.3 



Number of susy vacua with positive bosonic mass matrix 
Distribution of cosmological constants 

3.3.1 Density at zero cosmological constant 

3.3.2 Index density 

3.3.3 Total density 

3.3.4 Some general observations 

3.3.5 Distribution restricted to vacua with positive mass matrix 



3.4 Counting attractor points 



^ Distributions of nonsupersymmetric vacua 

[4.1| Supersymmetric and anti-supersymmetric branches 
Intermediate branches 
4.2.1 Large complex structure limit 
ITl Conifold limit 



4.2 



4.2.3 Distribution of cosmological constants in conifold limit 



[sj. Finding vacua with quantized flux 
|5.1| Flux vacua on rigid Calabi-Yau 







10 

HI 

13 

HI 

15 
16 
16 
17 
18 
19 

m 



29 

m 



36 

|38 



|6|. Conclus ions 



39 



- 1 - 



1. Introduction 



In this work, we study the distribution of metastable supersymmetric and nonsupersym- 
metric flux vacua in Calabi-Yau compactification of various string theories, along the lines 
developed in the works ||, 15] 

The effect of turning on gauge field strengths (or "flux" ) in string compactification has 
been studied in many works, starting with |^T[. Some recent examples include |2^, 25 1. 
Perhaps the most important qualitative effect of flux is that, since its contribution to the 
energy depends on the moduli of the compactification manifold, minimizing this energy 
will stabilize moduli, eliminating undesired massless fields. Since coupling constants in the 
low energy theory depend on moduli, finding the values at which moduli can be stabilized 
is an essential step in determining low energy predictions. It has also been suggested that 
taking into account the large number of possible choices for the flux, will lead to a large 
number of vacua with closely spaced values of the cosmological constant, and that some 
of these will reproduce its small observed value just on statistical grounds Q. Thus one 
would like to know the distribution of cosmological constants, and how this depends on 
the moduli and other parameters of the vacuum. 

There are a lot of flux vacua, and finding each one explicitly is a lot of work. Further- 
more, we are not entirely sure what properties we seek: there are many different scenarios 
for string phenomenology, each requiring different properties of the vacuum. Thus, rather 
than study individual vacua, we believe it is more interesting at this point to study the 
overall distribution of vacua in moduli space, and the distribution of quantities such as the 
cosmological constant and supersymmetry breaking scale. As discussed in |14, IC], such 
statistical results can serve as a guide to string phenomenology, and provide a "stringy" 
definition of naturalness. And, they are not much harder to get than results for individual 
vacua, as was seen in [|l|, |l5| and as we will demonstrate here. 

A useful way to state these problems is as that of finding vacua in a specific ensemble 
(or set) of AA = 1 effective supergravity theories, for which the Kahler potential and 
superpotential can be found explicitly. In principle, these theories are obtained by listing 
all string/M theory compactifications in a certain class, and in each case integrating out all 
but a finite number of fields, to obtain a Lagrangian valid at a low energy scale E. While 
not all vacua can be described by effective field theory, since at energies studied so far our 
universe seems to be described by effective field theory, this restriction seems adequate for 
the basic physics we want. 

One can certainly question whether this type of analysis captures all consistency con- 
ditions which vacua must satisfy. Perhaps the most important examples would be stability 
over cosmological time scales, and higher dimensional consistency conditions which are 
not obvious after integrating out fields. It is entirely possible that there are others; a list 
of speculations in this direction appears in Q, and these deserve study. Furthermore, it 
might turn out that early cosmology selects or favors a subset of preferred vacua. Our 
philosophy is not that we believe that none of this is important and thus can put absolute 
trust in the vacuum counting results below. Rather, we believe that, even if we had this 
additional information, it would not tell us which vacuum to consider a priori, and we 
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would still need to make an analysis of the type made here to find the relevant vacua, with 
the additional information taken into account as well. Thus the results we give should be 
considered as formal developments, with suggestive implications for real physical models, 
but which might be modified in light of better understanding. 

While the techniques we will describe could be used for any explicit ensemble of ef- 
fective supergravity theories, we work here with Calabi-Yau compactification in the large 
volume, weak coupling limit, because good techniques for explicitly computing and work- 
ing with the resulting supergravity Lagrangians exist at present only for this case. Indeed, 
type lib vacua on Calabi-Yau might turn out to be "representative" in the sense discussed 
in on various grounds. First, known dualities relate many other nonperturbative 

superpotentials to these cases, and it seems fair to say that all the structure in the potential 
which has been called on in model building so far, such as generation of exponentially small 
scales, and spontaneous supersymmetry breaking, can be seen in flux superpotentials. Sec- 
ond, dualities have been proposed which relate many of the other large classes of vacua to 
these. A systematic way to study the hypothesis that (say) lib on Calabi-Yau is represen- 
tative, would be to find statistics of vacua from two or more large sets of constructions; if 
both were representative, clearly these statistics would have to be the same. The present 
results are a necessary first step towards such a test, namely to find statistics for one large 
set of constructions. 

This concludes the justification of our approach. Our main discussion is somewhat 
technical, so we devote the remainder of the introduction to a basic overview of the type 
of results we will get. 

1.1 Distributions of vacua 

Our starting point is to imagine that we are given a list of effective supergravity theories 
Ti, T2, etc. all with the same configuration space (the space in which the chiral fields take 
values). We consider here theories with no gauge sector, so a theory is specified by a 
Kahler potential /Cj and superpotential Wi. 

We then apply the standard M = 1 supergravity formula for the potential, 

V = eW (^g^'WiWDjW* - ^1^1^) + ^^ (1-1) 

and look for solutions of dV/dz"^ = dV/dz^ = 0. Here Mp is the four dimensional Planck 
scale (which will shortly be set to 1). 

Vacua come in various types. First, as is familiar, a supersymmetric vacuum is a 
solution of 

DiW{z) = diW + ^{dilC)W = Q. (1.2) 

We consider both Minkowski = and AdS W ^ vacua. 

Conceptually, the simplest distribution we could consider is the "density of supersym- 
metric vacua," defined as 

dfisiz) =Y,Sz{DW^{z)) 

i 
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where Sz{f) is a delta function at / = 0, with a normahzation factor such that each solution 
of / = contributes unit weight in an integral J cP'^z. Thus, integrating this density over 
a region in configuration space, gives the number of vacua which stabilize the moduli in 
this region. A mathematically precise definition, and many explicit formulae which can be 
adapted to the physical situation, can be found in [15|. 
In this case, 

6,{DW{z)) = 6'^^''\DW{z))6'^'^\DW*{z))\detD^W{z)\ 

where the Jacobian term is introduced to cancel the one arising from the change of variables 
z — > DW. The matrix D^W is a 2n x 2n matrix 

= ( 9-.D,Wiz) d.D,Wiz)\ 

~ \ d-,DjW*{z) d^DjW*iz) J ' ^ ^ 

Essentially, this is the fermionic mass matrix. Note that we could have replaced the partial 
derivatives by covariant derivatives, as DiDjW = diDjW when DW = 0. 

With this definition, a vacuum with massless fermions counts as zero. There are better 



definitions, discussed in |14|, which would be appropriate if generic vacua had massless 
fermions. However, since generic vacua in our problem are isolated and have no massless 
fermions, this definition is fine. 

We can also define joint distributions such as the distribution of supersymmetric vacua 
with a given cosmological constant, 

dfis{z,A) = Y,Sz{DW,{z)) 6{A-{-3e'^^\Wi{z)\^)). 

i 

Below, we will define similar densities for nonsupersymmetric vacua of various types; at 
this point the basic idea should be clear. 

1.2 Approximate distributions of vacua 

Now, if we have a finite list of supergravity theories Tj, and if in each the number of vacua is 
finite, such a density will be a sum of delta functions. This is hard (though not impossible) 
to study, and for many purposes one might be satisfied with a continuous approximation 
to this density, a function p{z) whose integral over a region i?, 



(f'^zpiz), (1.4) 

R 



approximates the actual number of vacua in this region. 

What does this mean and what good is it? To give the question some context, suppose 
we had an explicit string theory construction of the Standard Model, and we were trying 
to decide whether it could reproduce the gauge and Yukawa couplings. While in some 
cases these are constrained by symmetry, this is not enough to determine the non-zero 
couplings. In the vast majority of explicit models, these couplings depend on moduli of 
the compactification (metric, bundle and brane moduli, etc.) and most results in this area 
address the problem of finding the formula for the couplings in terms of moduli. 
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Let us suppose we have such a formula. We could then use it to identify a region R 
in moduli space, or a region in the joint space of moduli, cosmological constant and other 
observables, such that any vacuum in this region would be guaranteed to produce couplings 
which agreed with observation to the required precision (let us say such vacua "work"). 
Thus the question would become, does the region R contain vacua, which can be obtained 
by stabilizing moduli. 

Now, suppose we had an approximation p{z, A) to the density of suitable vacua, with 
suitably small cosmological constant. Integrating it over the region R, would produce an 
approximation to the number of vacua which work. This starts to sound interesting, but 
of course we do not really want an approximate answer. In the end, we want to know if 
the approximation helps to answer the real question of finding a vacuum which works. 

If the region is small, or the density of vacua is small, one might need to interpret 
a result such as "approximately 10"^*^ vacua work." For the results we will discuss, this 
basically means that one expects that no vacua work, but it is possible that structure not 
reproduced by the approximation or some chance fine tuning will nevertheless lead to the 
existence of vacua. If there are competing classes of vacua which work, this would start to 
be evidence that vacua in the class under study do not work. To develop this hypothesis, 
the next question would be, what more do we need to do to prove that this region contains 
no vacua. 

If the region R is large enough to contain many vacua, there are two cases worth 
distinguishing. In our present state of ignorance, what we can typically do is try to enforce 
some but not all of the observational constraints on our model, and thus the number we 
would get at this stage would just be one factor in a final result. In this case, the most 
natural condition to put is that the integral Eq. (|1.4| ) should approximate the actual number 
of vacua N{R) with an error much less than N{R). 

Suppose we have solved the problem to the end, or we feel that the vacuum we have 
found is particularly interesting, say because it realizes some property of interest or refutes 
some conjecture in the literature. Then we would like to use our calculation of N{R) » 1 
to prove that a vacuum for the original, unapproximated problem, indeed sits in R; for a 
good approximation, this will be possible. 

Of course, we might decide that N{R) is so large, that the original goal of the discus- 
sion, say to test whether this class of vacua can reproduce the couplings in the Standard 
Model, becomes pointless. One can of course still hope that the assumptions that went into 
our choice of supergravity theories and definition of vacuum are false; however the compu- 
tation of N{R) under these assumptions would have been good enough and one would not 
need to improve that. 

This covers the various possibilities. The main points we want to make here are the 
following. First, what one wants to know next, and how one thinks about the problem, 
depends very much on whether N{R) >> 1, N{R) 1 or N{R) « 1, which is thus as 
important a question as finding particular vacua. 

Second, given that we trust our definition of "vacuum," it is actually less important 
to know how a given vacuum is obtained {e.g. by which choice of fiux) than to know that 



- 5 - 



it exists and stabilizes moduli (if the observables are controlled by the moduli, not by the 
fluxes; of course the fluxes could appear explicitly in the observables as well). 

Having said all this, we defer the discussion of how one could proceed in these various 
cases, to section 5 and to other work to appear. 



1.3 Other types of vacuum 

There are two types of nonsupersymmetric vacua we consider. The simpler possibility 
is breaking due to non-flux effects. In other words, we still seek solutions of Eq. (p..2|), 
calling upon other effects to break supersymmetry and lift the potential energy. This was 



invoked, for example, in [29|, which proposed to break supersymmetry by adding an anti 
D3-brane in lib compactification. Another possibility is to call on D term supersymmetry 
breaking, as has been discussed in many works. Indeed, to the extent the breaking can 
be understood in terms of effective M = 1 supergravity, this is the only possibility, and 
there are arguments in the literature that supersymmetry breaking by adding antibranes 
or misaligned branes is of this type [|^, ^, Thus, we are going to refer to this as D-type 
breaking. 

Before considering stability, the distribution of D-type vacua is the same as that for 
supersymmetric vacua, up to a factor which expresses the fraction of vacua which allow 
non-flux supersymmetry breaking. This might depend on the example at hand; we will 
simply set it to 1 here. 

Granting that in a given vacuum, adding the supersymmetry breaking term results in 
zero cosmological constant, we would identify the supersymmetry breaking scale as 

MLy = 3A 

where 

A = e'^("'^")|iy(z)|2 (1.5) 

is the norm of the superpotential for a vacuum stabilized at z. Of course, without the 
supersymmetry breaking, the cosmological constant of the resulting vacuum would have 
been AAdS = — 3A, so we will often refer to A as the "AdS cosmological constant." 

The other type of nonsupersymmetric vacuum is pure F type breaking; in other words 
to find a solution of y = which is not a solution of Eq. (|1.2| ) . The scale of this breaking 
is given by 

MLy = e'^g'^DiWD-W* 
which for y = is equal to the above. The density of these vacua is given by 



d^,.(z) = j;<5,(y'(z)). 



In general, one can have mixed D and F breaking. This is interesting only when both 
D and F terms depend on the same fields, which can only come about from non-flux effects, 
and is thus beyond our scope here. 
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In any case, the most interesting nonsuper symmetric vacua are the metastable (tachyon 
free) vacua, with V" positive definite. This constraint can be enforced formally by defini- 
tions such as 

dfiF,metastable{z) = ^ 5,{V' {z))e{V" (z)) 
i 

where 9{V") is 1 when the 2n x 2n real matrix of squared bosonic masses M = V" is 
positive definite. The derivatives of V appearing here are 

d^V = e^{DaDbWD''W - 2DaWW) (1.6) 
DadbV = e'^{DaDbD,WD^W - DaD^WW) (1.7) 

-2gb-aWW + DbD,WD-aD''W), (1.8) 

where R is the curvature of the cotangent bundle, i.e. R'^cab -^d = [Va, VjjXc = di{g^'^dagce) Xd- 
Note that we could have replaced the covariant derivatives by ordinary partial derivatives, 
because DdV = (fV when dV = 0. 

For D breaking, positivity M > can be analyzed as follows. First observe that in 
general, if DW = 0, 

M = H'^ - SA^/^F, (1.9) 

where 

H = 2d^k^''^. (1.10) 

This follows directly from Eq. ( |1.7] ) and Eq. (|l.§| ). Thus, to have M > 0, all eigenvalues A 
of H must satisfy A < or A > 3A^/^. In particular, if TV = at the critical point, M is 
automatically non-negative, and by continuity the same will be true for most susy vacua 
with small A. On the other hand, small positive eigenvalues of H will lead to tachyons and 
instability. 

The actual computations of all of these densities will of course rely heavily on specific 
details, but a general point worth keeping in mind is that any joint density 

dn{z,an) = '^6ziV'{z)) 5{an - Oniz)) 

i 

of vacua in moduli space along with any observables On{z) defined in terms of the Taylor 
series expansion of the effective Lagrangian about the vacuum (masses, couplings of moduli, 
etc.) can be computed if we simply know the joint distribution of W{z), IC{z,z) and a 
finite number of their derivatives evaluated at the point z, in other words a finite number 
of variables. Although obvious, this is very useful in structuring the problem, and is the 
main reason why this class of problem is so much simpler than problems involving flows 
on the moduli space. 

1.4 Index densities 

Finally, there is a quantity we call the "index density" . In the particular case of supersym- 
metric vacua, it is 

dls{z) = J2i-'^fUDWi{z)) 

i 



-7- 



where the index (—1)^ of a vacuum stabihzed at z is defined to be 



(-1)^ = sgn det D^W{z) (1.11) 

with D'^W{z) as defined in Eq. ([]|). 

The simplest reason to consider this is that the index is precisely the sign of the 
Jacobian which appeared in defining bz{DW^^ so 

{-\fb^{DW{z)) = 5(")(W(z))<5(")(W*(2))detL'2^(z) 

with no absolute value signs. Thus it is easier to compute, and provides a lower bound 
for the actual number of vacua. A formula for the index, and the explicit result for T^/Z2 
compactification, were given in [||]. 

There are also conceptual reasons to be interested in the supergravity index, as dis- 
cussed in ||T7| , 0]. To start, let us first comment on some differences between the problem 
of finding vacua in supergravity, and the much better studied problem of finding vacua in 
globally supersymmetric theories, satisfying 

diW{z) = 0. (1.12) 



Of course, Eq. (1.2) reduces to this upon taking the limit Mpi — > oo, or equivalently if all 
structure in W is on scales much less than Mpi (assuming derivatives of /C do not grow 
with Mpi). On the other hand, we need the supergravity correction to interpolate between 
different field theoretic limits. Of course, all hopes for getting a small cosmological constant 
out of Eq. (|l.lD rest on the supergravity term — 3|VKp as well. 



The problem of finding solutions of Eq. ( 1.12[) is holomorphic and therefore much easier 



than for Eq. (1^). In particular, vacua cannot be created or destroyed under variation of 
parameters, they can only move off to infinity or merge together. This makes it possible 
to give topological formulae for the total number of vacua in global supersymmetry; the 
possibility of vacua merging is accounted for by counting such vacua with multiplicity. 

Can we do the same for supergravity vacua? Evidently not, because one can construct 
a family of Kahler potentials JCt{z) such that varying t creates pairs of solutions of Eq. (|1.2|). 
One way to see this is to note that, in a region in which W ^ 0, Eq. ( |1.2| ) is equivalent to 
the condition that we are at a critical point of the function A from Eq. ( |1.5| ). Thus, where 
W ^ 0, one can apply Morse theory to this problem, as discussed in It is well known 
in this context that critical points can be created and destroyed in pairs. 

Clearly we cannot hope for a topological formula for the total number of vacua, but 
the above suggests using the Morse index for AAdS as a lower bound for the number of 
vacua, which might admit a topological formula. However, this is not correct because 
critical points of A^^^ are not necessarily vacua; indeed every point with W = (and /C 
nonsingular) is a critical point of A^ds- 

The search for a topological formula runs into other difficulties as well. Most impor- 
tantly, the configuration spaces which appear in known examples of effective supergravities 
are not compact, and cannot be compactified. The prototypical example is the upper half 
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plane. This configuration space has a boundary, the real axis, and it is easy to see in ex- 
amples that varying parameters (e.g. flux) can move vacua in and out of the configuration 
space. Furthermore, different flux sectors can contain different numbers of vacua. ^ 

Anyways, the correct generalization of the Morse index to this situation, as discussed 
in []T7| , |T^, is Eq. ( 1. 11 ). This agrees with the Morse index when W 0, and is (— 1)" 



when W = 0? 

bmce is a general (non-holomorphic) function, which away from W = can 

be deformed fairly arbitrarily by deforming /C, there is no obvious reason that one could 
not deform it to remove all cancelling pairs of vacua. ^ Thus, the index is the absolute 
minimal number of vacua which could be obtained by deforming the Kahler potential. 

Thus, one physical way to think about the index, is to say that the difference between 
the actual number of vacua, and the index, in some sense measures the number of "Kahler 
stabilized vacua," vacua whose existence depends on both the superpotential and the Kahler 



potential. We found in |15] and will find below that the actual density of vacua is typically 
the index density times a bounded function greater than one, so in this precise sense, there 
are many Kahler stabilized vacua. 

Since the DD terms in Eq. ( |1.3| ) go away in the limit Mpi oo, the index of a vacuum 
which survives this limit, and thus is not "Kahler stabilized," will necessarily be (—1)" (the 
same as for W = 0). Conversely, one could say that the vacua with index of the opposite 
sign are all Kahler stabilized, and would go away in this limit. One should realize however 
that this limit is highly ambiguous (the results change under Kahler- Weyl transformation) 
and that it may not in general make sense to say which particular vacua with index (— 1)" 
are Kahler stabilized or not. 

1.5 Ensembles of flux vacua 

We next discuss the set or ensemble of vacua we consider. General arguments have 
been given to the effect that the large volume, weak coupling limit of compactification 
of string/M theory on a Ricci-flat manifold M with fiux, can be described by a d = 4 
effective supergravity Lagrangian, whose configuration space C is the moduli space of com- 
pactifications on M with no fiux, the Kahler potential is taken to be the one for zero flux. 



and whose superpotential is the Gukov-Vafa-Witten superpotential [26|, which takes the 
form 



W= / GAn{z), (1.13) 

where z are the complex structure moduli, Q{z) is an appropriate form (depending on the 
theory and Af), and G is the p-form gauge field strength, which we normalize to have 
integral periods. 



^This statement is a bit imprecise; a more precise explanation taking duality into account is given in Q. 

^In |l| we instead used a convention for the index which includes an extra factor of ( — 1)", so that 
Minkowski vacua always count +1. There are arguments in favor of both conventions, and one should be 
careful to note which is in use. 

^The function IC must satisfy the constraint that ddIC is positive definite, but in one dimension this does 
not seem to prevent deforming away pairs of vacua. 
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Although we will discuss flux compactification of various theories: F theory on four- 
folds, heterotic string on CY3, and the (formally very similar) attractor description of black 
hole entropies, we work mostly with the lib flux compactifications on Calabi-Yau developed 



by Giddings, Kachru and Polchinski |22]. Then $7 is the holomorphic three- form on the 



CY, and G is a sum of the NS and RR three-form gauge field strengths 

with F,H e H^{M, Z), and r = C(o) + ie~'^ is the dilaton -axion. 

The "ensemble of effective field theories" with flux is then the set of supergravity 
theories with W given by Eq. ( |1.13| ), with G satisfying the tadpole constraint 

This is discussed in much more detail in 

We will leave out the Kahler moduli and forget about their contribution to the poten- 
tial, again for the reasons discussed in |jl|. We are going to derive many results relevant 
for this part of the problem, and will discuss it a bit in the conclusions, but reserve most 



of what we have to say about it to other work. |12]. 

With the above definition of W, the F-term potential V is given by [22|: 

V = 2Tse'^{DaWD''W -3WW), (1.14) 

where T3 = (27r)^^a'^^ is the D3-brane tension (in physical units). In subsequent sections 
we set 2T3 1, not as a choice of units (since Mpi = 1), but rather by shifting /C by a 
constant. 

Since Mp^ ^ = V^Mpi ^q, the dimensionless ratio V/Mpi^ ~ 1/^6^) but this is only 
because Mpi^4 grows with Ve. For orientation, in the traditional KK scenarios, Ve ~ ~ 
1/Mp; 10 up to 0(1) factors, so the natural energy scale of the flux potential is the string 
scale, and we will want supersymmetry breaking at scales Mgusyh < 1 or even << 1. In a 
"large extra dimensions" scenario, Vg/^f >> 1, and we might accept MsusyL ~ 1- 

Our basic results are obtained by neglecting the quantization of flux. This is expected 
to be a good approximation in the limit that the flux is large compared to other numbers 
such as the number of cycles. We will discuss this limit, and the sense in which the smooth 
distribution approximates the distributions of vacua at flnite L, in section 5. Our tools for 
doing this will be number theoretic theorems which state conditions on a "region in flux 
space" which guarantee that its volume provides a good estimate for the number of lattice 
points it contains. 

2. Notations and some useful formulas 

To avoid dragging along factors of in the calculations, we will slightly change notation 
in what follows and denote the usual holomorphic superpotential by W and reserve W for 
the Kahler invariant normalized but non-holomorphic superpotential: 

W{z,z) = e'^'^''''>/^W{z). (2.1) 
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Similarly we write 

n{z,z) =e'^^^^^^/^n{z) (2.2) 

for the normalized holomorphic form on the Calabi-Yau. We modify the definition of the 
covariant derivative accordingly: DaW = e^^'^DaW , etc. 

Consider first a general F-theory flux compactification on an elliptically fibered Calabi- 
Yau fourfold X. The Gukov-Vafa-Witten flux superpotential is 

W = [ G4An = N''Ua, (2.3) 



X 



where the IIq, = J A U are the periods of some basis {Sa} of H'^{X,'L). We normalize 
G4 such that G4 G H'^^X, Z),"^ so € 7}*. The flux has to satisfy the tadpole cancellation 
condition 

\ I G4AG4 = - iVD3 (2.4) 



27^"-'"- 24 

where Nd^ is the number of D3 branes minus the number of anti-D3 branes transversal to 
X. If we are looking for super symmetric vacua, this gives an upper bound 

L<L, (2.5) 

with L^: = xiX)/24:. Allowing anti-D3 branes, L^: can become bigger, but not indefinitely, 
as a sufficient number of anti-D3 branes in a flux background will decay into a state with 
flux and D3 branes only |27, 13 1. 



The superpotential depends on the complex structure moduli (a = 1, . . . , /i^'^(X)) 
only. The metric on complex structure moduli space is the Weil-Petersson metric, derived 
from the Kahler potential 

)C = -ln{Cl,n) = -In I nAh = -inUa{v~^T'^'ni3, (2.6) 
Jx 

where 7/„^ is the intersection form with respect to the basis {Sq}. 

We will often encounter intersection products of the covariant derivatives of the nor- 
malized period vector (or holomorphic 4- form): 

:Fa...b\c...d = Da--- Db^ r/-i Dc--- DdII = {Da--- Db^, Dc ■ ■ ■ Dd^ (2.7) 
^'a...b\c...d = Da--- DbUt]-^ Dc--- DdU = {Da--- Db^, Dc - - - Don), (2.8) 

where the capital indices can be either holomorphic or anti-holomorphic. These are most 
easily calculated by using identities of the form {DX,Y) = D{X,Y) — {X,DY), orthog- 
onality of (4,4 — A;) and (4 — /c',4)-forms with k' 7^ k, and commutation relations of D 
and D, together with Griffiths transversality, i.e. acting with k derivatives on the (4, 0)- 
form Q gives a sum of (4 — q, q)-ioTms with q at most equal to k. In fact, Da^ is pure 



*The flux quantization condition can actually be shifted by a nonintegral constant in some circumstances 
| p2| , but since we will make a continuum approximation for the fluxes anyway, we can ignore such subtleties. 
Also, in F-theory, not all fluxes in _H"'(X, Z) are allowed: essentially, one leg should be on the elliptic fiber 
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^ T.\h — 


-9ab 


•r- 

ah\cd 


^acbd ' 


T , 

•' ab\cd 


9ba9cd 


T' — 

•' ab\cd 




Jx 



(3,1) and DaDi^O, is pure (2,2), which can be shown in similar fashion. As an exam- 

pie, we have T,^-, = {Dan,D-,n) = Da{n,D-,n) - {n,DaD-,n) = - {n,g^j^n) = -g^i. 

Similarly, the absence of a (3,l)-part in DaDb^l is follows from J^ab\c = {DaDf,Q, Dc^) = 
Da{Di,Q,, Dc^) — {DifQ, DaDc^) = —Dagbc~ {Db^,gac^) = 0. Thus, most lower order 
J^-tensors vanish. Some nonzero ones are 

(2.9) 

^ 9bd9ac + Qaddbc (2.10) 

(2.11) 

nAdadbdcdan =J^abcd- (2.12) 

It should be noted that in general DaDbD^fl is not pure (1,3); there can be a (2,2)-part, 
since {DaDbDc^, D^De^) = DaJ^bc\de — ^aRbdce-' which in general is nonvanishing. 

2.1 Orientifold limit 

Things simplify considerably in orientifold limits of the F-theory compactification. Then 
X = (T^ X y)/Z2 with Y a Calabi-Yau threefold, which is equivalent to type lib on the 
corresponding orientifold of Y with constant dilaton-axion r. 

2 1 

There are n = {Y) complex structure moduli of Y surviving the orientifold projec- 
tion, and 21^2^ {Y) + 2 fluxes can be turned on |0, Let = t and (i = 1, . . . , n) be 
the complex structure moduli of Y . Then 

= fii(t°) AJ]3(i*), (2.13) 

so /C = /Ci + /C3 with /Ci = — ln(z(Oi, ^i)) = — ln(2 Im r) and /C3 = — \n{i{Q^, ^3))) and 
the metric and curvature components mixing and i all vanish. In general, there may 
be other fourfold moduli as well, which take X away from the orientifold limit. These 
correspond to D7-brane moduli from the lib point of view. We will ignore them in what 
follows. 

As before, we define normalized holomorphic forms by ^lr = eF"'l'^tlr. Using the same 
methods as we used before to compute intersection products of derivatives of 1^4, one 
obtains 

DqVLx = J'o^i, DoDoQi = (2.14) 
with Tq = DqQi) = ie'^^{fl,i,dofti) = — 1/(t — r), and 

DiD^Qs = J'^jkD^^i, (2.15) 

with J^ijk = DiDjDk^'s) = ie'^'-^ {(l^^didjdk^^) ■ Therefore all components of J^ab\cd 
and J^abcd are zero, except 

^Oi\Oj = 9oo9i~j (2.16) 

^ ij\kl = ^ij^^mkl (2-17) 
^Qijk = ^O^ijk- (2-18) 



^The minus sign refers to the part of the cohomology odd under the orientifold involution. 
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Note also that DiDj^l^ = J^Qij^D^D^n^. 

The space of allowed fluxes Hp(X) consists of harmonic 4-forms G4 = —aAF^ + fdAH^, 
where {a, /?} is a canonical basis of harmonic 1-forms on (such that ili = 13 — ra) and F3 
and are harmonic 3-forms on Y (identified with type lib R-R resp. NS-NS flux). The 
main simpliflcation occurs because 17, Da^l, DoDiO, and their complex conjugates form 
a Hodge-decomposition basis of Hp{X). To see this, note that we have dmiHF{X) = 
2(2n + 2), which equals the number of vectors in the proposed basis set, and that linear 
independence of this set follows from the intersection products computed earlier. The basis 
can be turned into an orthonormal basis by introducing an orthonormal frame e\ for the 
metric on moduli space, (5^^ = e^^^ge^, where capital letters refer to the frame indices.^ 
We take Cg = J^q, so J-q = 1. The basis B = {Q., Da^, DqDjCI} U {c.c.} now satisfies 

(e,^) = diag(l,-l„+i,l„,l,-l„+i,l„). (2.19) 

Various physical quantities have a simple expression in terms of components with respect 
to this basis. Writing 

G4 = Xn-Y^DA^ + Z^DoDiQ + cc, (2.20) 



we get for example for the flux superpotential Eq. ( |2.3| ) and its derivatives (transformed 
to the orthonormal frame by Da " " " = • • • Da ■ ■ ■ D),): 

W = {G4,n)=X (2.21) 

DaW = {G^,DAn) = YA (2.22) 

DoDgW = (2.23) 

DoDiW = Zi (2.24) 

DjDjW = TukZ"" (2.25) 

D^_DjDjW = TijkY"^ (2.26) 

DiDjDkW = {DiTjkl)Z'^ + TijkY\ (2.27) 

for the potential 

y = |y|2 _3|x|2, (2.28) 
and for the flux induced D3-charge tadpole 

L = \^Xr]N = ^(G4, G4) = |Xp - |y|2 + \Zf. (2.29) 
3. Distributions of supersymmetric vacua 

A supersymmetric flux vacuum is characterized by a choice of K flux quanta A^" and a 
solution to DW = 0. We wish to compute the total number of such flux vacua satisfying 



'For explicit numerical indices 0, 1, ... we will underline frame indices, but only if confusion could arise. 
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the constraint L = ^Nr]N < L^,, 

Nsusy{L<L,)= ^{L-L,) (3.1) 

susy vac 

= J_ / ^e-W(a) (3.2) 

where C runs along the imaginary axis passing zero to the right, and where we introduced 
the Laplace transformed "weighted number" of vacua 

AA(a)^ J^e-t^^^ (3.3) 

vac 

= V/ (f'^z5'^'^{DW)\detD'^W\ (3.4) 
AT Jm 

^[ (f'^z [ d^Ne-"^^"^^ 6^'^{DW)\detD^W\. (3.5) 

Jm J 

In the last step we approximated the sum over fluxes by an integral. By rescaling — > 
N/y/a^ it is easy to see that Eq. (|3.5D scales simply as a~^/^, so in this approximation 
Eq. dU) gives: 

L 

Nsusy{L < L,) = e{L,)j^^N{a = 1). (3.6) 

As discussed in the previous section, in the orientifold limit we have m = n + 1 with n 
the number of complex structure moduli of Y, and K = 4m. In this case, it is possible to 
directly evaluate the Gaussian integral by changing variables from to (X, Y, Z, X ,Y , Z), 
related to each other by the Hodge decomposition Eq. (|2.20|): 



AT = r]-^{XU - Y^DaU + Z^DqDjU + c.c). (3.7) 

The Jacobian for this change of variables is 

J = 22'"|detM| = 4"*|detr/|-^/Vet(Mtr/M)|^/2^ (3.8) 

where M = ?7~^(n, —D^II, Do-C/II, c.c). The extra factor 2^™ accounts for the fact that 
for complex variables we use the convention (Pz = ^dz A dz. Happily, because of the 
orthonormality of our Hodge decomposition basis B, we have 

M^r]M = diag(l, -l„+i, 1„, 1, 1„), (3.9) 

hence the Jacobian is simply J = 4"^! detr/|~^/^. Furthermore, from Eq. (2.22), we get 

6'^"'{DaW) = \dete^\-^d^"'{DAW) = {det g)-^ S^"" (Ya) (3.10) 

and from Eq. (|2^ -Eq. (|]25|) together with DaDiW = g^^W: 

{detg)-^detD^W = det[{Do, Dj, Dq, Dif ■ {DqW , DjW , DqW, DjW)] (3.11) 



-14- 



det 



\ 



det 



ixrdet 



( X Zj \ 
bijX Zi TukZ"" 
Zj X 
\ZiTijkZ'' bijX ) 

/ X Zj 
X Zj 
Zj 6ijX :FjjkZ^ 
\Zj TukZ"" 6jjX J 

( bjjX - ^ TjjkZ"^ 

\ TjjkZ^ 5jjX-^ 



(3.12) 



(3.13) 



(3.14) 



X 



Putting everything together, we find for the total number of supersymmetric vacua: 

(27rL,)2™ 



M{L < L,) 



(2m)! 



det7?ri/2 / d^'^z det gp{z) 



(3.15) 



M 



where 



TT 



~2m 



d^Xd^^Ze-l^l'-l^I'lXPldet 



ZtZ 



( hjX 

\ TjjkZ^ bjjX 



X 
K 



TjjkZ'' 



(3.16) 



X 



The function p measures the density of supersymmetric vacua per unit volume in moduli 
space. It is specified entirely in terms of the special geometry data J-jjk- In particular, 
p has no dependence on the dilaton modulus r, and therefore the integration over the 
fundamental r-domain in Eq. ( 3.151 ) simply contributes a factor vo1(7V1t-) = 7r/12. 

Similarly to p, we define the index density Pmd-, counting vacua with signs, by dropping 
the absolute value signs from the determinant in Eq. (|3.16| ). 

3.1 Computing densities 
3.1.1 The case n = 1 

The total susy vacuum number density for n = 1 can be computed explicitly from Eq. ( |3.16| ): 

P = ^'^l d'Xd^Ze-l^l'-l^l' \\X\^ + \Z\' - (2 + |.F|2)|X|2|Z|2| 

3^ - (2 + |J^|2)rs| 



TT 



drds e ^ * + 



(3.17) 
(3.18) 
(3.19) 



This can be obtained by splitting up the integration domain in three parts, separated by 
the lines s/r = ^(2 + \ J-\'^ it IJ'^I y^4 + |J^p) on which the determinant changes sign. The 
first two terms in this expression correspond to the index density: 



Pind = (2 - |.F|')/vr2 



(3.20) 



In the large complex structure limit one has universally \J-\ = 2/\/3 and, in the spe- 
cial coordinate t, = — 3/(t — t)^ (to verify the former, recall that !F = {e\YTttt = 
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{gtt)~^^'^e'^Hnr]-^dfn, with e'^s = -^^^ and Ur]-^dfU = 6k). So pics = S/vr^ (and 
Pind = yo/3), and if we approximate the large complex structure region M.LCS by the 
standard fundamental domain in the upper half plane, we get 

d^rgrf [ dhgup « -L x ^ x 2 = i-. (3.21) 

Mr JMlcs 

Near a conifold point (or more generally a discriminant locus), T blows up (see example 
2 below). Note that in that case, up to a sign, the total number density equals the index 
density. In fact this is true for any m: if all Tijk oo, the terms involving Tijk in 



Eq. (3.16) will dominate the determinant, so 

deiD'^W ^ {-iy'\X\^\deiJ^ijKZ^\^ (3.22) 
and putting the absolute value signs around D^W only removes the overall (—1)"". 
3.1.2 Index density 

Computing the total density for n > 1 becomes hard. However, the index density can 
be given a simple expression in terms of geometric quantities [Q]. One way to do this is 
to rewrite the determinant as a Gaussian over Grassmann variables, and then to per- 
form first the Gaussian over X and Z, and next the the Grassmann integral. Using 
^IJKL = ^IK^^MJL - Sij5kl - ^iL^jK (which follows e.g. from comparing Eq. (|2.1C| ) 
with Eq. ( 2.17] )) and i?gooo = — 2, one gets after some manipulations 



dpind = d'"^zdetgpind = vr"'" det(i? + u;l), (3.23) 

with R the curvature form and uj the Kahler form on Ai. This is in agreement with earlier 
results [|l|. Note however that the index density does not factorize, det{R + u;l)j'2y.Y 
det(i? + uj1)'j^2 a det(i? + uj1)y, so the claim in vl of ||l[ that adding the dilaton to the 
moduli just multiplies the index by 1/12 was not correct. For example for n = 1, one has 

dM^nd = 7r-2detf^° + ^° + '^^ ° ] (3.24) 

= TT-'^i-LOo+uJi) A{Ri+uJo+uJi) (3.25) 
= -vr^^wo Ai?i. (3.26) 

In this case, we therefore simply have 

dpind = -^xiMY). (3.27) 
M -"-^ 

3.1.3 Example 1: 

As a toy example, let us take Y to be the {T'^)^/Z,2 orientifold with the and the 
fluxes restricted to be diagonal and symmetric under permuations of the three T^'s. Then 
the complex structure moduli space is the fundamental domain in the upper half plane, 
parametrizing the modulus, and Eq. ( |3.21| ) is exact. The orientifold has 64 03-planes, 
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so = 16. A basis for the symmetric fluxes is {'^a}a=i,-,'i, with the given by the 
generating function '^'^=1 ^at°'~^ = Y\k=ii'^k + tPk)- Here {ak,Pk) is a canonical basis of 
H^{Z) of the A:th T^. The nonvanishing intersection products on are (Si,E4) = 1 and 
(S2, S3) = 3. To avoid a subtlety with flux quantization involving discrete fluxes on the 03 
planes |20|, we will furthermore as in l2|] restrict to even fluxes, i.e. we take as basis {2Sq, 
The corresponding intersection form on Y thus has as nonzero entries r/^ = — = 2 and 
r/l^ = — ?y32 = 6. The intersection form on x y is the direct product of this with the 
intersection form eij. Therefore | detry| = (2^ x 9)^, and the total number of these flux 
vacua is, according to Eq. ( 3.15[) : 



(27r X 16)'' , 4 , 1 1 , , 

Msusy = ^j-^ X (2^ X 9)-i X — = 1231. (3.28) 



3.1.4 Example 2: Conifold 

Let y be a Calabi-Yau manifold near a generic conifold degeneration. For simplicity we 
only consider one modulus, namely the period of the vanishing cycle v = Jj^^- Inclusion of 
more moduli does not change the essential features. The monodromy around v = implies 
that the period of the dual cycle is of the form J^Cl = — ^ v + analytic terms. The metric 
near f = is then 

g^^ ^clnf- (3.29) 

\v\ 

where /i is some constant and c = e^°/27r with /Co the Kahler potential at w = 0. Further- 
more 

/ \ / 2 \ -3/2 

= g-^l'' e'^ iij^ Cl 53 J^Cl + anal. J « J c In J ^ , (3.30) 

so as announced earlier, we see that J- ^ 00 when v ^ 0. The same is true for p ~ |.Fp/7r^. 
However, the density integrated over the fundamental r-domain and < i? remains finite. 
For small R: 

cfrg-rf I cPyg^apM ^-3-. (3.31) 

121n^ 



Note that the constant c has dropped out of this expression. Plugging this in Eq. ( 3.15| ), 
we get for the number of susy vacua with L < and \v\ < R: 

^vac = *-T- (3.32) 

181n^ 

The logarithmic dependence on R implies that a substantial fraction of vacua are extremely 
close to the conifold point. For example when L,, = 100 and ^ = 1, there are still about 
one million susy vacua with |f | < 10"^'^''. Interestingly, vacua very close to conifold degen- 
erations are precisely the desired ones in the context of phenomenological model building, 



as they provide a natural mechanism for generating large scale hierarchies [22|, and may 



enable controlled constructions of de Sitter vacua by adding anti-D3 branes, as proposed 



by KKLT [29|. However, for the latter it is also necessary that the mass matrix at the 
critical point is positive, and as we will see below, this condition dramatically reduces the 
number of candidate vacua. 
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3.1.5 Example 3: Mirror Quintic 

The mirror quintic is given by a quotient of the hypersurface + + X3 + + = 
5^p X1X2X3X4X5 in CP^. It has one complex structure modulus, ^p, whose fundamental 
domain A^y is the wedge — 7r/5 < argV' < vr/5. Its periods are well known Q and can be 
expressed as Meijer G-functions, which makes it possible to study this case numerically. 



Ttpg /12 




-3-2-1 12 3V 



Figure 1: The susy vacuum number density per unit ip coordinate volume, npg^^ /12, on the real 
■f/j-axis, for the mirror quintic. 

Fig. H shows a plot of ^ pg^^, i-e. the susy vacuum number density per unit ip coordi- 
nate volume, on the real ip-axis (the factor 7r/12 comes from integrating over r). The drop 
for I'^l > 1 is due to a similar drop in g^^; p itself tends to the large complex structure 
value when tp ^ 00. The divergence at i/' = 1 is due to the presence of a conifold 
singularity there. In the notation of example 2, the parameters specifying g and near 
the conifold are p = 8.94 and c = 1.26 x 10"^ (with v « -Air'^b-^/'^ {ip - 1)). 

We numerically computed^ the integrated susy vacuum number density. We found: 

dp = 5.46 X 10"^ (3.33) 

M 

This can be compared with an estimate of the large complex structure contribution, ob- 



tained similar to Eq. (3.21) by using the LCS expressions for g and p, but now cutting off 
the integral say at Im t = 2. (Here t defined by bip = e"^'^**/^, and the conifold point is 
located at Im = 27r/51n5 = 1.28.) The result is l/167r = 1.99 x 10"^. The exact 

numerical result for this region is almost the same: 1.97 x 10^^. Thus for orientifolds of 
the mirror quintic, about 36% of all susy flux vacua are at Im t > 2 (and this fraction is 

^This was done as follows. First, we divided the moduli space in patches, since different regions have 
different suitable coordinates and special care is required near singularities. In each patch, the periods and 
their derivatives were evaluated on a dense grid of points, and an approximation of these functions was 
constructed by interpolation (because direct evaluation of Meijer functions is much too time-consuming). 
Finally from this data the various desired quantities were constructed and integrated. 
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proportional to one over the lower bound on Im t). On the other hand, using Eq. (|3.31j ), 
we get that the fraction of vacua with 1-0 — 11 < 5 <C 1 equals 0.486/ ln(6. 41/5^). For 
S = 10-3, this is about 3 %, and for S = lO-^o still 1%. 
For the integrated index density, we found 



dfii^a = -1.666 X 10"^ PS -1/60 (3.34) 

Combined with Eq. ( p. 27 ), this indicates that x(-^y) = 1/5- Indeed, this can be verified 
analytically. The integral of the Euler class can be written as a sum of boundary contour 
integrals as follows: 



27r X 5 



9dlng^^ (3.35) 
My 

Bdlng^^ (3.36) 



where the Pi are the 5 copies of the conifold point in the V'-plane. For if; ^ oo we have 



9^^- = 3/4|^|^ln^|V'|, 



so 



ain9** = -i(l + i^)* (3.38) 

and the corresponding contour integral produces a contribution (z/lOvr) x (— 27ri) = 1/5 to 
Eq. (|0^ ). On the other hand, near the conifold point tp — 1, g^^ — cln(//V|V'-lp), and 

dlng^^p = -- ^2 T— T^^' (3.39) 

In V - i 

so the corresponding contour integral is zero. Adding up all contributions, we thus see that 
X = 1/5. 

Numerical integration of the volume of A4 gives 5vol(A4) = 3.1416 « vr. Again, this 
can be understood topologically, using 

uj = ^ddlC = ^ddlni{n,n), (3.40) 

plus the fact that if Cl is normalized such that /C is regular at ip = 0, we have ~ 
I/IV'P In^ \ ip\ for ip ^ oo. Writing the volume integral as a sum over contours, again only 
the ip = CO contour contributes, and this contribution equals vr, as expected. 

3.1.6 signature distribution 

For some applications, such as the RG flow interpretation of domain wall flows, one needs 
to know whether the critical point of is a maximum, a minimum, or a saddle point of 
\W\. This corresponds to a positive, negative or indefinite Hessian = {da,daY • 
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(c^lM^I, SftlVFl). At a critical point, one has 9a9fe|T^| = ^i^DaDbW and dadb\W 



1 w 

2 \W\- 



DaD^W, so from Eq. ( 3.12| ) we see that we have to investigate the eigenvalues of 



d^\W\ 



21X1 



/|X|2 


xz 

\ xz TXZ 



XZ \ 
xz TXZ 

|xp 
|xp / 



(3.41) 



Clearly the sum of the eigenvalues tr(i^|VF| > 0, so there are no maxima (this is true in 
general, since 5a5f;|T^| ~ and g > 0). In general, a matrix is positive definite iff 

all upper left submatrices have positive determinant. In the case at hand, this implies the 
conditions + |Z|^ - (2 + \T\^)\X\^\Z\^ > and |Xp - jZp > 0. This restricts the 
integration domain of Eq. ( |3.18D to one of its three segments, and thus the density of susy 
vacua which are minima of |iy| is 



(3.42) 



More information can be obtained by looking directly at the eigenvalues. We just quote 
the results: = and 



1 



(3.43) 



----- ^2 yjT^- 

In the large complex structure limit one has, rather democratically, = = 

= p/3. In the conifold limit on the other hand, = p and = = 

0. 

3.2 Number of susy vacua with positive bosonic mass matrix 

Due to the properties of AdS, supersymmetric vacua are always perturbatively stable, even 
if the critical point of the potential V is not a minimum. Obviously, this is no longer true 
if supersymmetry is broken and the cosmological constant is lifted to a positive value. In 
particular, if as in the KKLT scenario supersymmetry is broken by adding an anti-D3 brane 
to a supersymmetric AdS vacuum, thus shifting the potential up by a constant^ such that 
the critical point of the potential becomes positive, the original AdS critical point should 
be a minimum in order for the lifted vacuum to be perturbatively stable. It is therefore 
important to compute the number of supersymmetric vacua which are local minima of V. 

A critical point is a minimum if the mass matrix d^V = {d, BY ■ {dV, dV) is positive 
definite. Using Eq. (|]|)-Eq. (|l|) and Eq. (^)-Eq. (|^, we get, at a supersymmetric 
critical point (i.e. Y = 0): 



DodoV = 
DidoV = -ZiX 



(3.44) 
(3.45) 



^constant in the complex structure deformation directions 
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DodjV = -ZjX (3.46) 

DidjV = -TijkZ'^X (3.47) 

D^doV = -2|X|2 + (3.48) 

Djd^V = TiklZ^'z'^ (3.49) 

D^djV = TjklZ'^Z'^ (3.50) 

DidjV = -2|X|2 + ZiZj + .F7i^*^.FMjLZ^Z^ (3.51) 

We can use covariant derivatives of V here instead of ordinary derivatives because dV = 
implies DdV = d'^V . 

Let us work out the case n = 1, which has 



/|Z|2-2|X|2 jrz2 -XZ \ 

(1 + |jr|2)|^|2 _ 2|X|2 -XZ -J^XZ 
-XZ |Z|2-2|X|2 Jc-22 

\ -XZ -/"XZ /■Z2 (1 + |J^|2)|Z|2 - 2|X|2y 

(3.52) 

A matrix is positive definite iff all its upper left subdeterminants are positive. With r = 
|X|2, s = |Z|2, this gives the following conditions: 

s - 2r > (3.53) 
4r2-2 (2 + |J^P) rs + s^ > (3.54) 
(s-2r) (4r2 - (5 + 2|J^|2) rs + s^) > (3.55) 
(16 - 4 (2 + I r s + s^) (r^ _ (2 + | J^|2) rs + s^) > 0. (3.56) 

A straightforward analysis of these inequalities shows that they boil down to simply 



|Z|V|X|2 = s/r > 4 + 2\T\^ + 2|J^| ^4 + |J^|2. (3.57) 

To compute the number density of susy vacua with positive mass matrix, we should there- 
fore evaluate the integral Eq. ( |3.18| ) with (r, s) restricted to the region satisfying this 
condition. This gives 



+ 179 + 42 \J^\^ + 96 y^I+JW + 42 l-Fp ^/I+JW ,o .on 

PM>o = ; (3-58) 

5 + 2|^|2 + 2|.F| y/I+W) 7r2 

as shown in fig. ^. At J'^ = 0, the relative fraction of susy vacua with positive mass matrix 
is approximately 39% and maximal. In the large complex structure limit = 2/\/3 the 
fraction is about 19%, and in the conifold limit it is zero. Indeed, while the total density p 
grows quadratically with \ J^\ when ^ oo, the density pm>o decreases quadratically with 
\!F\. This means that very near a conifold point, though there are many susy vacua, only 
an extremely small fraction has positive mass matrix! 

Let us make this more precise. In the notation of our example 2 in the previous section, 
we have for small v. 



PM>o^^(ln^y\v\\ (3.59) 
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2 4 6 8 |F| 

Figure 2: The upper curve shows the total susy vacuum density as a function of |JF|, the lower 
curve the density of susy vacua with positive mass matrix. 



Integrated over r and |f | < R, this gives 

.2 



J cfrgrf J d\ g,^ pm>o = |^(3 + 67 + 67' + 47^ + 2-f^)R\ (3.60) 



,,2 



where 7(-R) = In^. Because of the R dependence, this rapidly goes to zero with R. 
If we take the parameters of the mirror quintic conifold for example, i.e. /i = 8.94 and 
c = 1.26 X 10-2, and we take L^, = 100, then the expected number of susy vacua with 
M > drops below 1 for \v\ < 0.02, or j-^ - 1| < 0.004. Increasing to 1000, these 
numbers become just one order of magnitude smaller. Thus, for n = 1 and with reasonable 
parameter values, at most only a modest hierarchy of scales can be generated through the 



conifold throat mechanism of |22], if we insist on having a positive mass matrix. 

On the other hand, if a near-conifold vacuum has positive mass matrix, Eq. ( 3.57| ) 
together with \X\'^ + |Zp = L < shows that it automatically has a small value for 
\W\ = \X\, and therefore a small cosmological constant. We will make this more precise in 
the next section. 

The positivity properties oi M = d?V for susy vacua can also be analyzed as follows. 
First observe that in general, if DW = 0, 

M = H'^ -3\W\H, (3.61) 

where 

H = 2(f\W\. (3.62) 



This follows directly from Eq. (|l.7| )-Eq. (|l.8| ). Thus, to have M > 0, all eigenvalues A of 
H must satisfy A < or A > 3\W\. In particular, if Vl^ = at the critical point, M is 
automatically non-negative, and by continuity the same will be true for most susy vacua 
with small W. The suppression of vacua with M > near a conifold point (for n = 1) can 
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now be seen in the following way. According to Eq. ( |3.41| ), the matrix H can be written 
as H = \X\1 + AH, with 



= U n ' S = — {^^-] (3.63) 





The eigenvalues of AH are (±Ai,±A2), with AfA^ = det = \Z\^ and 2(Af + A^) = 
tr (AH)'^ = 2(|.Fp + 2)|Zp. When — > oo, the eigenvalues of AH are therefore approxi- 
mately given by ib|J'^|^"'^|Z|, and the eigenvalues of -ff by A = \X\ ±\J^\'^^\Z\. The condition 
on A to have M > translates to \Z\ > 2\J^\\X\, in agreement with what we found earlier. 
Note that the "dangerous" eigenvectors (the eigenvectors of AH with eigenvalues ~ 
are approximately aligned with the dilaton direction, i.e. {1,0[1/J^],1,0[1/J^]). In 
a sense, the special form of the matrix Eq. ( |3.63| ) leads to a sort of "seesaw" mixing with 
the dilaton, and the small eigenvalue. This may be specific to n = 1; a similar analysis for 
n > 1 suugests that there is no longer suppression of M > vacua near generic points of 
the discriminant locus. 
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Figure 3: The number density of susy vacua with positive mass matrix, per unit coordinate volume, 
on the real i/^-axis, for the mirror quintic. Compare to fig. |l|. 



As an example, a plot of the number density of M > vacua per unit coordinate 
volume is shown in fig. |^ for the mirror quintic, on the real -i/^-axis. The sharp dip near 
'0 = 1 is due to the conifold singularity. For the integrated density we find 

/ diiM>Q = 5.12 X 10-^ (3.64) 
JM 

Comparing this to Eq. ( |3.33|) , we thus see that about 9% of all susy vacua has positive 
mass matrix. This is not too far from the naive guess 1/16 based on the fact that there 
are four mass eigenvalues, with each a 50% chance to be positive. 
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3.3 Distribution of cosmological constants 

The cosmological constant in a supersymmetric vacuum is = — 3|M^p (in 2T3 = 1 units). 

We wish to count the number of susy vacua with L < and 1^ > K = — 3A*, i.e. with 
< A*. Analogous to Eq. (^^, we have 

MsusyiL < L,, \W\^ < A,) = [ cP'^z f^*dX ^ [ — e"^*i/(z, a, A) (3.65) 

Jm Jo 27ri J a 

where 

i^{z,a,X) = J d^Ne-''^6{\W\'^ - \)6^"\DW)\detD^W\ (3.66) 

= a^-^''^ j d^N e-^ 5{\W\'^ - aX)5'^'^{DW)\deiD'^W\. (3.67) 
Parallel to Eq. ( glsl) , we can write Eq. ( |3.65| ) also in the form 

Msusy{L < L,,\W\'' < A,) = ^^^^^|detr/ri/2 / d^^z^ [^*dXpiX,z) (3.68) 

(-^/^j! Jm Jo 

The quantity dXp{X, z)/p{z) then gives the fraction of vacua at z with \ W\'^ in a width dA 
interval around A. 

Specializing again to the orientifold limit, so ii' = 4m = 4(n + 1), we make the change 
of variables from to {X, Y, Z). In these variables the additional constraint in Eq. ( |3.67| ) is 
simply |Ap = a A, which we can solve by putting X = \/aX, since the integral is invariant 
under an overall phase transformation of {X,Y,Z). Thus we find, similar to Eq. ( 3.1(^ ) 
(but with the additional constraint): 



(.,a,A) = -^p^/d-Ze---l-ndet 



/ VaX Zj \ 
SijV^ Zj TukZ^ 
Zj \faX 
V Zi TijkZ^ 6ijV^ J 



2m- 



/m 
d^~^Ze-\^\\Y^Cu{,Z,Z){aXf\ (3.69) 

k=0 



where g = detg, rj = |detr/| and the Ck{Z,Z) are homogeneous polynomial functions 
obtained by expanding the determinant. For example in the n = 1 case, 

^.(^, a, A) = e-"^ [ d^Z e-\^\" \{aXf - (2 + \T\'')\Z\^aX + \Z\\ (3.70) 

^a-^ J 

3.3.1 Density at zero cosmological constant 

For physical applications, the most interesting quantity is the number of vacua at very 
small cosmological constant, i.e. the case A=k <C -L*. To a good approximation, we can then 



^We could also have implemented the < A» inequality by a similar Laplace transform, but here it 

is slightly easier to solve the constraints directly. 
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neglect the higher order terms in Eq. ( |3.69| ) and simply compute the density at A = 0. The 
integral becomes 

I(^)./.-Ze-I^P|det(°^^j^;_,.)|=. (3.71) 

This is a significant simplification compared to Eq. ( |3.69| ), as this can be evaluated using 
Wick's theorem or by rewriting the determinant as a Grassmann integral and then doing 
the integral over Z, similar to what was done to obtain Eq. ( |3.23| ). Doing the integral over 
a in Eq. ( p.65| ) now gives 

27rm 



P\\=o 



vr 



(3.72) 



Note that for n = 1, this is independent of J- and therefore of z: po = S/vr^L*. Applying 
this to the example of the mirror quintic, for which vo1(A^t- x M^) = f > we get for 
small A*: 



Ar{L<L.,\Wf<X.) = ^^ 

45 



(3.73) 



For = 100, this becomes 10 ''A*, so the expected smallest possible cosmological con- 
stant is |A| ~ 10~^T3. For L^, = 1000, this is three orders of magnitude smaller. 




Figure 4: Left: Vacuum number densities (true and absolute value of index) in the large complex 
structure limit {!F — 2/V3), as a function of cosmological constant value. Right: Same near the 
conifold limit (.F = 100 for this example). Note that despite appearances, p is non-zero at A = as 
expressed in Eq. ( ^.73 ); it just does not receive the conifold enhancement there. 



3.3.2 Index density 

The distribution for arbitrary cosmological constants is harder to compute, mainly because 
of the absolute value signs in Eq. (3.69). Let us therefore drop these for now, so we count 



vacua with signs. Then after integration over Z, we will still have a polynomial in a: 



Vindiz, a, A) 



k-2m+l 



(3.74) 



fc=0 
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and therefore, by comparing Eq. (|3^ ) and Eq. 



In the case n = 1, we thus get for the cosmological constant density counted with signs 

P^nd{z, A) = -^(1 - x)(2 - (10 + 3|^|2) X + (14 + 3|^|2) 0(1 - x) (3.77) 

where we wrote x = As noted above, the density at A = equals S/L^.tt'^. In the large 
complex structure limit J- 2/-v/3, Pind = 8(1 — x)(l — 7x + 9a;^)/L*7r^ and in the conifold 
limit F ^ CO, for x not too small, Pmd/I-^P = —12(1 — x)^x/L*7r^. 

3.3.3 Total density 



Let us now turn to the true density. The absolute value signs in Eq. ( 3.69 ) make the 
integral hard to evaluate in general, but for n = 1 one gets something of the form 

u{z, a. A) = Po(aA)e-"^ + P_(aA)e-''-"^ + P+(aA)e-^+"^ (3.78) 

with the Pi polynomials and 6_, 6+ some ^-dependent coefficients (see below). Doing the 
a-contour integral, this then gives, explicitly: 

L.vr^ p(z, A) = L,^Vmd(2, \)+P{x){e{l-b_x)-e{l-h+x))+Q{x){e{l-b_x)+e{l-h+x)) 

(3.79) 

where 

P{x) = -4 + 6(4 + x - 12(4 + x^ + (32 - 6|.F|^ - \T\^) x^ (3.80) 

Q(x) = (4 + |.F|2)3/2i73^3 (3_g^) 



h± = -74Tm^(v/4Tm^± |.F|) (3.82) 
In the large complex structure limit 

16 

L^tt'^p{z, A) = L^Tr'^pindiz, A) + 16(l-4x)2 (9(l-4x) + — (16x-3)(3-4x)2 6'(3-4x), (3.83) 

and in the conifold limit, for x not too small, 

p{z, A)/|.F|2 ^ 12(1 - xfx/L,TT^ = -p,nd{\ z). (3.84) 

This is not immediately obvious from the above expressions, but can be seen fromEq. ( |3?70| ): 
the middle term dominates, and the absolute value just removes the minus sign. As dis- 
cussed before, the convergence of index density and true density can be expected to hold 
more generally near degenerations where J-ijk — *■ oo (see also below). 

These considerations are illustrated for the large complex structure and conifold limits 
in fig. ^. Note that the approximation of the true density by the index density is perfect 
near the conifold, and near the extremities of A, and qualitatively still not bad away from 
those. 
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3.3.4 Some general observations 

A few simple observations can be made about the general case. It is obvious that the 
cosmological constant always satisfies < -L*. This is because \W\'^ + |Zp = L < L^: hy 
construction. The density at A = is zero for the same reason; in fact, the index density 
near that point is suppressed by (L — X^:)"^~^, as follows from Eq. ( 3.75| ). By contrast, the 



density at A = is always nonvanishing. Note however that in case all J-ijk oo, this 
value is much smaller than the total density p{z) (compare Eq. ( 3.7l| ) to Eq. ( [3.221) ). In 



general, for large m, the distribution can be expected to peak at a small value of A, because 
of the suppression factor {L — \^)'^~^ . More intuitively, the lower \W\'^ is, the more "room" 
there is for Z to satisfy \ W\'^ + \Z\^ < L^. 

Finally, in the limit where J^ijk — > oo, we can compute p{z^ A) more explicitly, since 
then, for A not too small, 

" (Sl^2^ / 17^"*^ / ^^"^e-"^-'"'^«A|det(^,..Z-)P (3.85) 
= Axil-xf"^-^ (3.86) 

where x = A/L,, and A = ^^^^ / f^^"^ e"'^'' | dei{J^ijKZ^)? . This distribution has 
mean fj, = (x) = 2/ {2m + 1) and standard deviation = ((x — /i)^) = (2m — l)/(m + 
l)(2m + 1)^, so for large m, most susy vacua in this limit have a cosmological constant A 
of order —T^L^/m. This should be compared to the string scale energy density, which in 
Einstein frame is Tp = gs{2iTa')~'^ = 2TTgsT^- Therefore in this limit vacua with |A| well 
below the string scale are comparatively rare, unless m S> L^/gg. 

3.3.5 Distribution restricted to vacua with positive mass matrix 

If we restrict to vacua with positive mass matrix, Eq. ( |3.57| ) together with |Zp + |Ap = 
L < L*, implies the cutoff A < L,/(5 + 2|J^|2 + 2|J^|y^4 + |J^|2). In the large complex 
structure limit, this is A < L*/13, and near the conifold limit, A < L*/4|.Fp. Therefore, 
positive M vacua near the conifold point automatically have small cosmological constant 
(small compared to the string scale) . Recall however that positive M vacua are suppressed 
near the conifold point. We saw earlier that for the parameter values of the mirror quintic 
and = 100, no positive M vacua are expected below v 10^^. At this point ^ 20 
and thus A < 10^'^L,^ = 0.1, hence the cosmological constant cutoff is not much below the 
string scale. Increasing to 1000 decreases the A cutoff with just one order of magnitude. 
Thus, for n = 1, the positive M vacua closest to the conifold point are not expected to be 
close enough to force the cosmological constant to be hierarchically smaller than the string 
scale. 

This does not mean of course that there are no M > vacua with very small cosmo- 
logical constant; in fact, the vacuum density at A* = is the same for M > 0- vacua as 
for vacua without constraints on M, since if X = 0, the M > condition Eq. ( p. 57 ) is 



automatically satisfied. Most vacua with very small A will therefore have a positive mass 



matrix. This is actually true for any supergravity theory, as follows from Eq. (1.9). 
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3.4 Counting attractor points 

The techniques we used to count supersymmetric flux vacua can also be used to count 
super symmetric black holes, or attractor points, in type IIB theory compactified on a 
Calabi-Yau X. More precisely, we wish to count the number of duality inequivalent, regular, 
spherically symmetric, BPS black holes with entropy S less than S^,. The charge of a black 
hole is given by an element Q = A^"Sq, of H'^{X, Z). The central charge Z = {Q, $73) plays 
a role similar to the (normalized) superpotential W: the moduli at the horizon are always 
fixed at a critical point of \Z\, i.e. DZ = {d + dlC)Z = 0. The entropy is then given by 
S = A/4: = tt\Z\'^, evaluated at the critical point. For a given Q, critical points may or 
may not exist, may or may not be located in the fundamental domain A4, and may or may 
not be unique, but clearly we can label all equivalence classes of black holes by a charge 
vector N and a critical point within the fundamental domain. Thus our desired number 
is, in a continuum approximation (valid for large 5"*) similar to what we had before: 

AfsniS < S,) = — I ^e"^*/" f d^N [ d^^ze""!^!' ^^"(DZ) | det D^^j. (3.87) 
2Tn J a J 

Here K is the number of flux components and n the number of complex structure moduli, 
K = 2n + 2. After rescaling N N/y/a, doing the a integral becomes straightforward, 
and we get 

AfBHiS < S,) = ^j^^ [ d^^zdetgpiz) (3.88) 
(n + 1)! 

where 

Piz) = j^^^y^^T^g J d^-^^Ne-\^n^-iDZ)\detD^Z\. (3.89) 



To evaluate this integral, we change variables to a Hodge-decomposition basis of H'^{X), 



similar to Eq. ( p.2[)| ): 

Q = iXn - iY^Djn + C.C., (3.90) 

so Z = {Q,0,) = X and DjZ = {Q,DiZ) = Yj. Capital indices again refer to an 
orthonormal frame. The Jacobian for the change of variables from to {X,Y,X,Y) can 
be computed in a way analogous to what we did for flux vacua; the result is J = 2""*"^|r/|~^/^, 
where \r]\ is the determinant of the intersection form (equal to 1 if we sum over the full 
charge lattice). Furthermore, using the special geometry identity DjDjZ = TjjkY^ , we 
get that at a critical point (where Y = Q): 

\D,D.,Z DiDjZ ) \ a SijXj ' ' ^ ' 

Plugging this in Eq. ( |3.89| ) gives 

p = [ d^X [ d^T e-l^l' 6^'^{Y) lAl^'^ (3.92) 

„M/i J J 



n|„|l/2 ■ 



7r"|r/ 



(3.93) 
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The factor det g dropped out because of the change from coordinate to orthonormal frame. 
Note that this expression for p is independent of z. This means that attractor points are 
uniformly distributed over moduh space. Our final result is, for large 5*: 

2"+ivol(A^) 

For the mirror quintic, n = 1 and vol(A^) = vr/5, hence 

AfBHiS < S,) = ^5,2 (3.95) 

The problem of counting the number of duality inequivalent black holes with given 
entropy was studied in |30| using number theory techniques. In particular, for Y = T^xK3, 



this was shown to be related to class numbers and the number of inequivalent embeddings 
of a given two-dimensional lattice into the charge lattice. 

To compare to our results, we restrict to algebraic K3's, as the complex structure 
moduli moduli space of a generic K3 is not Hausdorff, so the volume factor in Eq. ( 3.94| ) 
would not make any sense. An algebraic K3 has moduli space 

MK3 = 0{A)\0{2,k)/0{2) xO{k), 

where A is the charge lattice, which is the orthogonal complement of the Picard lattice and 
has signature {2,k). This space is Hausdorff and has complex dimension k. A black hole 
charge is specified by two charge vectors p,q £ A. 

In this case, counting black holes with given entropy, or more precisely with given 
"discriminant" D = (S/n)'^ (which is always an integer for K3 x T^), amounts roughly 
to computing the number of inequivalent lattice embeddings in A of lattices spanned by q 
and p with determinant D}^ This number can in principle be obtained from the Smith- 
Minkowski "mass formula" , which was used in [^0| to derive an estimate for the asymptotic 
growth of the number N(D) of inequivalent black holes with discriminant D. The result 
is N{D) ~ d'^/^ and consequently the number of black holes with D < will grow as 



To compare with Eq. ( 3.94| ), note that n = A; + 1, so we get 



Mbh{D <D,) = ^-^^o\{Mt2)^o\{Mkz) dI-^^'^, (3.96) 

which agrees with the growth given above, but is a bit more precise. Turning things around, 
this formula should give a predicition for the asymptotic behavior of the Smith-Minkowski 
mass formula. 

4. Distributions of nonsupersymmetric vacua 

A general flux vacuum, supersymmetric or not, is characterized by a flux vector and a 
critical point of the corresponding potential V , which for now we do not require to be a 



"We thank Greg Moore for explaining this to us. 
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local minimum. For nonsupersymmetric vacua, the condition L < L* is no longer sufficient 
to guarantee a finite volume of allowed fluxes. Physically, it is reasonable to put also a 
bound on the supersymmetry breaking parameter: = \Y\ < F*. Such a bound, 

together with the bound on L = |Xp — |yp + \Z\'^, is indeed precisely what we need to 
make the volume of allowed fluxes finite. 

In our approximation, the total number of flux vacua satisfying these bounds is then 
given by 

J\f{L<L^,\DW\ <F^) = 2^'^ j d^^z j (fX(f'^(f''Ze{L^-L)e{F^-\Y\)6'^'''{dV)\deicPV\ 

(4.1) 

where (fV = (5,9)* • {dV,dV), and using Eq. (H)-Eq. (|r|) and Eq. (|2^ ~Eq. (|^: 



doV = ZiY^ - 2YaX (4.2 

diV = ZiY^ + J^ukY-^Z^^ - 2YiX (4.3 

DodoV = (4.4 

DidoV = TiklY^Y^ - ZiX (4.5 

D^djV = TjklY'^Y'^ - ZjX (4.6 

DidjV = 2TijkY'^Y'' - TukZ'^X + DiTjklZ'^Y'^ (4.7 

DqBqV = -2|X|2 - 2|yo|^ + iXfP + \Zi^ (4.8 

Di%V = -Yi% + TiklZ'^Z'^ (4.9 

D^djV = -YoYj + ^jklZ^'z'^ (4.10 
DiBjV = -2|X|2 + |yoP - 2YiYj + ZjZj + Tik^' Tmjl{Y^Y^ + Z^ Z^) (4.11 

We can use covariant derivatives of V in the integral instead of ordinary derivatives because 
if dV = 0, then DdV = d^V . 

4.1 Supersymmetric and anti-supersymmetric branches 

The main complication arises because the constraint dV = is quadratic. It defines a 
cone in C™, which has several branches. One obvious branch is Ya = 0, with X and 
Zj arbitrary. This corresponds to the supersymmetric vacua discussed in the previous 
sections. Another obvious one is X = Z/ = 0. As we will see, the vacua on this branch 
behave in a way as "anti-supersymmetric" vacua. For example, while susy vacua have 
imaginary self-dual fluxes, these have imaginary anti-self-dual fluxes. There are more 
complicated "intermediate" branches too, which arise when the constraints considered as 
linear equations in {X, Z) (or in Y) are degenerate. 

On the supersymmetric branch, we have 5'^'^{dV) = | det M|~^5^™(y), with M = 
{dy^dyY ■ {dV,dV). This threatens to make the integral divergent when detM 0, 
but note that by the chain rule and because Y = d^V\Y=o = {d,dY ■ {dV,dV) = 
{d, BY • (y, y) • (%, Oy)* • {BV, dV) = D^W-M, so the factor det M cancels out of Eq. (|3) 
and we are left with the integral we had before for the supersymmetric case, as we should 
of course. 
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Similarly, on the branch X = Z = 0, we have S^'^idV) = \det M\-^6^'^{X, Z), with 
M = {dx,dz,dx,dzY ■ {dV,dV). Again detM cancels out, because (fV\x=z=o = A- M, 
with 



A = {Do,Dj,Do,Djy ■ {W,DoDjW,W,DoDjW) 

/ Yj Yo \ 
6ijYo Yi J^ijkY^ 
Yo Yj 
\YiJ^ijkY^ 6ijYo J 



(4.12) 
(4.13) 



so the remaining determinant factor in Eq. (4J.) is 



det^ = llol^det 



'6ijYo 



iXj. 



51) 

J'ijkY'^ 



TijkY"^ 
^uYo-^ 



(4.14) 



which is identical to the det D^W factor Eq. (3.14) of the super symmetric case, after 
substitution X ^ Yq, Zj ^ Yj. If we moreover replace the supersymmetric constraint 
L = |Xp + < by |yp < (or, since now L = -\Y\'^, by L > -L*), the 
integral for this nonsupersymmetric branch is therefore formally the same as the one for 
the supersymmetric branch! Thus 



x=z=o^ 



\DW\ < F,) = MsusyiL < Fl 



(4.15) 



In fact this can be understood directly from the defining equations. If (F^, H^jT, z^) with 
Im T > solves F^ — tH^ = i * (F3 — tH^), which is the condition for a supersymmetric 



vacuum, then {F^, H'^,t' , /') = (-F3, R-^ 



-r, z 



has Im r' > and solves F!^ - 



-t'H' 



-I* 



(F3 — t'H^), which is the condition for a nonsupersymmetric vacuum on the branch under 
consideration. An equivalent way to see this map is to observe that Dq^Nrr + tN^s) ~ 
—]^RR — T^NS together with the fact that the action of Dq and Dq interchanges {X, Z) 
and Y. Note however that this is not a map between vacua with the same topological data, 
since it maps L —L; the susy vacua have positive and the nonsusy vacua negative L. 

The cosmological constant of these nonsupersymmetric vacua is A = 2T3|L'Typ = 
— 2Lr3. Since L is quantized, this is always at least of the order of the string scale energy 
density Tp = IngsT^, so actually the field theory approximation on which this analysis is 
based cannot be trusted, and the existence of these vacua in the full theory is doubtful. 



4.2 Intermediate branches 

The anti-supersymmetric branch X = Z = the constraint cone dV = is parametrized 
by the values of the F-terms Y^. For generic values of the Y^, the unique solution to 
dV = is indeed X = Z = 0, since it is just a generic, complete system of linear equations 
in (X, Z). However, for some values of Y the linear system can become degenerate, namely 
when detM = where M is as before given by M = (5x, 5^, 9jf , 9^)* • {dV,dV). This 
happens at the intersection with other branches. 
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We will only analyze the case n = 1 here. Then 

doV = ZiYi - 2YoX (4.16) 
diV = ZiYo +TY1Z1 - 2YiX (4.17) 

and detM = -4(|Fo|'' + - (2 + |.?='n|loPl^iP)- The equation detM = has two 
branches of solutions: 

|YiP = Xl\Yo\^ with A± = = ± V4H^). (4.18) 

These branches can be parametrized for example by (Yq, Zi). The explicit solutions of the 
constraints dV = are then, for generic {Yq, Zi): 

Yi = A± e*('^'-g-P'-2^'-s^i) Yo (4.19) 

X = ix^eii'^''sF-'i>^rgZr+2>xrgYo) (4 20) 

Solving the constraint in this way, the delta-function d'^{dV) produces a factor IdetMl"-*^ 
with M = {dY^-idx-idYii^xY ' [dVidV). On both branches we have 

|detM| =4|jr|V4+|jr|2|yo|'|^i|'. (4.21) 

Furthermore, in general, 

<iV^(ff). (4.22, 

where for n = 1 

- 2|X|2 - 2|yoP + l^iP ^ZI-YqYi \ 



R 



TZl - loll l^ll' - 2|X|2 + |FoP - 2|yip + |J^|2(|yi|2 + |Zi|2) 



/ ty:1-xz^ \ 

\ty:1- xzx t^tyoY^ -xzi) + d^tY^Zi j ' 

4.2.1 Large complex structure limit 

The general expression for de\,(PV is obviously very complicated, but at large complex 
structure, where J- = 2/\/3 and D!F = 0, things simplify considerably. In this case, we get 
for the integration measure on the first branch, with r = |lbP) s = \Z\\^: 

IdetdVl/ldetMl = ^|(12r - 5s)(8r - 3s)| (4.23) 

and on the second branch 

25 

|det(iV|/|detM| = — |(24r - 7s)(4r + 3s)|. (4.24) 
972 

The flux number L, susy breaking parameter |y| and cosmological constant V for the first 
branch is (still in the large complex structure limit) 

L = -4r + ^ (4.25) 

|yp = 4r (4.26) 

F = 4r - ^ (4.27) 
4 
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and for the second one 



4r 13s , 
L = -- + — (4.28) 



It is most convenient to implement the constraint on L in Eq. (4J) by doing a Laplace 
transform, and the one on \ Y\'^ by constraining r directly. The result for arbitrary F^: and 
L* is a bit complicated and not very instructive — it involves several polynomials multiplied 
by different step-functions — but for positive and not too large, the step functions 
are unimportant, and we get a simple polynomial. On branch 1, for x = / < 6: 

Mi^icsiL < L^,\DW\ <F,) = x (3.07- 0.499 x + 0.0256x2 + 1.04 x^)yLt (4.31) 

where V = vol{Mr x Mies) = (7r/12)vol(A^/cs). On branch 2, for x < 42/25: 

J^2,ics{L < L*, \DW\ <F^)=x (5.59 + 0.416 x - 2.31 x^ + 1.88 x^) V L^. (4.32) 

This can be compared to the number of supersymmetric flux vacua at large complex struc- 
ture, 

^fsusy{L<L,) = ^^^^^^V = 13.2 V Li (4.33) 



We still have to investigate if the vacua are perturbatively stable. For vacua with positive 
cosmological constant, this is the case if and only if (PV is positive definite. A matrix 
is positive iff all upper left subdeterminants are positive, which here translates in the 
conditions 

2r-s>0 (4.34) 
-(12r-7s)(4r-3s) > (4.35) 
-(2r-s)(96r2-68rs + 15s2) > (4.36) 
rs(12r-5s)(8r-3s) > (4.37) 



on branch 1 and similar conditions on branch 2. Straightforward analysis of these systems 
of inequalities shows that they have no solutions, on either branch. Alternatively, we could 
have computed the eigenvalues of d?V. For the first branch for example, these are 

{4^,4^_^,_^,_8r + 3s} (4.38) 

Again we see it is impossible to have only positive eigenvalues. Thus we find, somewhat 
surprisingly, that there are no meta-stable de Sitter vacua at large complex structure on 
these branches. 
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4.2.2 Conifold limit 

In the conifold limit, we have, as in Eq. ( 3.30| ): 

i 



(4.39) 



where 7 = In -j^. Now DT is no longer zero: 



i(l - 3/7) i(7 - 3) 



07^ ip' 



(4.40) 



However DT turns out to drop out of all relevant quantities to leading order in e. 

Prom Eq. ( |4.19| ), it follows that on branch 1, to leading order, |Yi| = |lo|/|e| and 
\X\ = \Zi\/2\e\, and similarly but with opposite power of e on branch 2. On branch 1, we 
thus have, to leading order: 



L 



\Y\ 



(|Zi|V4-|loP)/|6| 

|2 



klV4-|y|' 



l^olVkP 



\y\ 



V 



\y\ 



3|z|V4 



(4.41) 
(4.42) 
(4.43) 



Since the constraints in Eq. ( [4.1D keep L and |y| finite, Zi and Yq are of order e, and 
the rescaled variables y,z are of order 1. In these variables the integration measure is, to 
leading order in e: 

det d?V 



detM 



,2 ,2 \y?\z? 
d yd z . 
16 er 



(4.44) 



Doing the integral over y and z (implementing the constraints as before) yields for Eq. ( [4.1| ) 



= ^ [LiQ{U) - [U - 3F2)(L, + Fl)H{U + Fl^ I d'Td\^. (4.45) 



Apart from the L^,-dependent bracket, this is the same as the expression for the number of 
near-conifold super symmetric vacua, so using Eq. ( |3.32| ) we immediately get for the number 
of nonsusy critical points on branch 1 with \v\ < R, L < L^, and \Y\ < F^,: 



vr 



181n^ 



{Lt9{L,) - (L, - 3F2)(L, + F^f9{L, + F^)) 



For positive, the L,,-dependent bracket reduces to QL^F^ + 8L^F^ + 3F^. 
The same computation can be done for branch 2. To leading order in e: 



The final result is identical, A/i 



L 
V 

det (fV 



l^il'-|>ol' 



\y\ 



3|e|2|z|V4 



detM 

on — 2,con 



\y\ 

,2 ,2 \y\W 

a yd z — |— — . 



(4.46) 



(4.47) 
(4.48) 
(4.49) 

(4.50) 
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Do these near-conifold nonsusy critical points have positive mass matrix? On the first 
branch, the first subdet condition is 2r — s > (again with r = |yp and s = and 
using this, the third condition becomes 

(16 + 0[e^])r^ + (2/|e|2 + 0[e°])r s + {I + 0[e^]) < 0. (4.51) 

For smah e, this is obviously false, so (fV can never be positive definite on this branch 
in the near-conifold region; Mi^con,M>o = 0. For the second branch on the other hand, 
the positivity conditions are, after some simplifications and dropping terms which are 
negligible^^ for all r, s at sufficiently small e: 

s/r > 4/|e|2 (4.52) 
-72s2 + 4rs/|e|2+7siner-^/2^^/2|g| ^ ^4_53) 



where 6 = arg.F-|- argYi — argZi and 7 as under Eq. ( 4.39 ). With s/r = 4ii^/|ep, this 
becomes 

u-1 > (4.54) 

p(n) = 7sin6lu^ - 7^ + 1 > 0. (4.55) 

The polynomial p{u) has three real roots Uj. We have p'{0) = 0, p"(0) = —27^ < 0, 
p{0) = 1 > and p{l) < 0. Therefore, if sin 6 < 0, then tii < ^2 < < ^3 < 1 and p{u) < 
for u > 1, so the system of inequalities has no solutions. When sin 9 < on the other hand, 
we have ui < < U2 < 1 < U3, so the above system of inequalities boils down to 

u > U3 ^ ^/ sinO. (4.56) 

The latter approximation becomes exact in the extreme conifold limit 7 — > 00, but is 
already very good for 7 > 4 (error less than 0.5 %). This inequality together with the 
constraint L < implies |yp = r < ^ , so the super symmetry-breaking parameter 



F = \Y\ is automatically less than e-/Ll/2^. If we take the susy breaking cutoff greater 
than this number, the integral Eq. ( |4.1| ) becomes independent of -F* and is given by 



2 .con.Af >0 



j StcPv deigv{v), (4.57) 



^ e sm 

167r " 



v(v) = ^f-^ I de 1'' dr I dsrs (4.58) 
2 kr ./n ./n . 



(4.59) 



12874 

Up to logarithmic factors, this has the same dependence on v as the number of supersym- 



metric near-conifold vacua with M > 0, which according to Eq. ( |3.58| ) is given by 

l^susyiv) = (4.60) 

o 

Nonsusy near-conifold vacua with M > are sparser than their susy counterparts; their 

density ratio is 1/147^ <C 1. 

^^This includes putting 7 + 0[1] ~ 7. It is possible to keep terms of lower order in 7, but this only 
complicates the formulas without changing the essential features. 
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4.2.3 Distribution of cosmological constants in conifold limit 

The distribution of cosmological constants for M > vacua near the conifold point is 
obtained by adding 6{V — A) to the integrand of Eq. ( 4.58| ) , with V as given by Eq. ( 4.4g| ) : 



i^{v,A) = ^ [ de [ dr [ * dsrs6{r-^s-A) (4.61) 

2|er Jo Jo J^PX- ^ 

sin^ 

2 riT rM 

de / drr{r-A) (4.62) 



IGvr / 4 \ 

V3?y 



JO 



where M = min(-^^^^^, + A, The first two entries in min(. . .) come from 

the integration boundaries of s. Since s > 0, M has to positive, and we get the following 
condition on A to get a nonzero result: 

_ ^ < A < 0. (4.63) 

In particular we see that the cosmological constant can only be negative for these vacua, 
so there are no meta-stable de Sitter vacua in the conifold region either. This is because 
the M > condition forces the positive term in V to be smaller than the negative one. 
Under these conditions on A, we have in fact that M = — ^3™ ^ , except in a very small 
interval (width of order 1/7^ relative to Eq. ( [4.63| )) near the lower bound on A. Outside 
that interval, the integral is therefore straightforward to compute and equal to 



Inside the small boundary interval, the density quickly drops to zero. 



5. Finding vacua with quantized flux 

So far, we have been discussing the problem of finding the volume of the region in "flux 
space" which contains vacua, not the problem of counting physical vacua with quantized 
flux. This approximation was necessary at a very early stage, when we solved the conditions 
DW{z) = or V'{z) = for some of the fluxes in terms of the others, as at a generic point 
z in moduli space these conditions will have no solution with integer fluxes. 

Nevertheless, we can get results for the problem with quantized fluxes using these 
techniques. The approach will be to consider a region R in moduli space, and character- 
ize the corresponding region Sl in the "space of fluxes" which contains vacua with 
N'qN/2 < L. Since the number of flux vacua in this region of moduli space is the number 
of points with integral coordinates, i.e. points in Sl n Z^, we need ways to estimate this 
number from facts about the geometry of Sl. 

The region S containing supersymmetric vacua is not hard to describe . At a given 
point (z, r) in moduh space, the 2n + 2 = K/2 conditions = iV^ReAII" = iV^ImAII" 
pick out a subspace Sz,r — ^^^"^ of the space of fluxes. The region S is then the union of 
the individual Sz,t for all values of moduli {z, r) S R. One can show that the subspaces 
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vary transversally (the "non-degeneracy condition" of Q) and thus this wih fih out a K- 
dimensional region. Finahy, the subset 5^ C 5 is obtained by imposing the additional 
constraint NrjN < /2L. 

It is clear that S" is a cone; in other words AS* = S. 
Furthermore, for a small region R, the constraints DW = 
will not vary much, so one can think of S as roughly a 
cone over the product of a K/2 — 1-dimensional sphere (at 
fixed radius in M^/^) with a i^/2-dimensional ball. Over 
a small region, the constraint Nr]N/2 < L is not very 
different from a positive definite quadratic constraint, so 
Sl is roughly the r < \fL part of this cone.^'^ For example, 
for K = 2^ one can have regions Sl as pictured in figure 

i 

Having characterized Sl, our goal is to count how 
many lattice points it contains. The most basic estimate is 
of course that, if we make L sufficiently large, this number 
of lattice points will be well approximated by the volume 




of Sl- Furthermore, it is intuitively clear (see [21|, 2.iv 
for more precise statements) that the leading correction 
to this is proportional to the surface area of the boundary 
of Sl, so for large L 

N{L) = L^/V(5i) + L^'/2-l/2^(^^) + 0(^i^/2-l/2-.) 

(5.1) 

where V{Si) and A{Si) are the volume and surface area 
at L = 1. 

How big need L be to reach this scaling regime? In addition to the condition L^/'^V{Si) » 
1 for the region to contain many lattice points, one clearly needs 

A(5i) 



Figure 5: Two possible regions 
S and S' containing supersym- 
metric vacua in a if = 2 flux 
space. S contains many fewer 
quantized flux points at small L 
than its volume, while S' contains 
many more. 



VI » 



ViSi 



(5.2) 



as well. 

To illustrate what happens for smaller values of L, consider the case of K = 2, and the 
two cones S and S' illustrated in figure ^. Cone S, misaligned with the lattice, does not 
contain any lattice points near the origin. On the other hand, cone S' , aligned with the 
lattice, contains roughly r = \/L (positive) lattice points within a small distance r from 
the origin, far more than its volume V ~ r'^9/2 (where 6 is the opening angle). 

This phenomenon persists all the way out to ~ 1, at which point the two dimensional 
nature of the cones starts to be visible, and the estimate Eq. ( |5.lD becomes valid. This 
leads to the condition r > 1/6. Since V{Si) = 6/2 while ^(^i) ~ 2, this is the same as 
Eq. (P). 

Now the quantity V{Si) is just the integrated vacuum density over the region R, and 
the quantitv A(S^ ) can be comput ed as an integral over moduli space using the same 
^^Sl is however not convex, and the theorems of Minkowski and Mordell called upon in vl of this paper 
are not apphcable. 
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techniques. In the approximation that the vacuum density is just the volume form on 
moduh space, the surface area of the boundary will just be the surface area of the boundary 
in moduli space. Taking the region to be a sphere in moduli space of radius r, we find 



so the condition Eq. (|5.2|) becomes 



L > 



K 



(5.3) 



Thus, if we consider a large enough region, or the entire moduli space in order to find 
the total number of vacua, the condition for the asymptotic vacuum counting formulas we 
have discussed in this work to hold is L > cK with some order one coefficient. But if we 
subdivide the region into subregions which do not satisfy Eq. (^]^), we will find that the 
number of vacua in each subregion will show oscillations around this central scaling. In 
fact, most regions will contain a smaller number of vacua (like S above), while a few should 
have anomalously large numbers (like S' above), averaging out to Eq. ( ^.1[ ). 

5.1 Flux vacua on rigid Calabi-Yau 

As an illustration of this, consider the following toy 
problem with i^T = 4, studied in ||]. The configuration 
space is simply the fundamental region of the upper 
half plane, parameterized by r. The flux superpoten- 
tials are taken to be 

W = At^B 

with A = oi + ia2 and B = h\^ ih^ each taking values 
in Z + iZ. This would be obtained if we considered flux 
vacua on a rigid Calabi-Yau, with no complex structure 
moduli, 63 = 2, and the periods Hi = 1 and 112 = i- 
The tadpole condition N'qN/2 < L becomes 



One then has 



lmA*B < L 



DW = ^ f 



B 



(5.4) 



(5.5) 




-0 . 5 



. 5 



Figure 6: Values of t for rigid CY 
flux vacua with Lmax = 150. 



Thus, it is very easy to find all the vacua and the value 
of T at which they are stabilized in this problem. We 
first enumerate all choices of A and B satisfying the 
bound Eq. (|5.4[) , taking one representative of each orbit 
of the 5L(2,Z) duality group. As discussed in [||], this can be done by taking 02 = 0, 
< 61 < oi and 0162 ^ L. Then, for each choice of flux, we take the value of r from 
Eq. (^]^) and map it into the fundamental region by an SL{2, Z) transformation. The 
resulting plot for L = 150 is shown in figure 0. 
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A striking feature of the figure is the presence of holes around points such as t = ni 
with n € Ij. At the center of such holes, there is moreover a big degeneracy of vacua. For 
example there are 240 vacua at r = 2i. This clearly illustrates the phenomena discussed 
above. Starting with a tiny disk D around r = 2i, the corresponding cone in flux space 
is very narrow, but aligned with the lattice so it captures 240 points. For a somewhat 
displaced disk, the cone is not aligned with the lattice and captures no points. Increasing 
the radius of the disk makes the cone wider, and at a certain radius new lattice points 
enter. 

Despite the intricate structure of the finite L result, it is true that a disc of sufficiently 
large area A will contain approximately 2'kAL'^ vacua. This is illustrated for L = 150 
in figure 0, where estimated and real numbers of vacua are compared in discs around 
r = 2z of stepwise increasing radius. Note that the first additional vacua enter the circle 
at a coordinate radius R ~ 0.12, and that just beyond that radius, at i? ~ 0.15, the 
approximation becomes good. The corresponding radius in the proper metric is r ~ 0.04. 
Since the holes around the integers clearly correspond to the worst case scenario for the 
estimate, we can thus conclude empirically that for L = 150, our estimate will be good 
when r > 0.04. Moreover, by comparing the results for different values of L, we found that 
the radii of the holes scale precisely as 1/^/L. Thus, our empirical result for the reliability 
of the approximation is 



r > 



0.5 

7z' 



which is compatible with Eq. (|5.3|) . 
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Figure 7: Number of vacua in a circle of coordinate radius R around r = 2i, with R increasing in 
steps dR = 0.01. Pink bars give the estimated value, green bars the actual value. 



6. Conclusions 

We gave general results for the distribution of supersymmetric and nonsuper symmetric flux 
vacua in type lib string theory, and studied examples with one complex structure modulus 
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in detail. Let us conclude with a brief summary of the results, some comparisons to what 
one might expect intuitively, and questions for further work. 

A simplified picture of the results is that one can define an "average density of vacua" 
in the moduli space, which can be integrated over a region of interest and then multiplied 
by a "total number of allowed choices of flux," to estimate the total number of vacua which 
stabilize moduli in that region. This estimate becomes exact in the limit of large flux, 



and should be good for flux satisfying the bounds discussed in section 5, say Eq. (|5.3D for 
supersymmetric vacua. 

The zeroth approximation for the average density of vacua is the volume form on 
moduli space constructed from the metric which would appear in the effective supergravity 
kinetic term; explicitly if this is 



then we find 



in F theory and lib flux compactiflcation, with K B4 fluxes and tadpole bound L. 

Of course, the actual results we obtained were of course more complicated, depending 
not just on the metric but on curvatures and its derivatives. However it is still useful 
to think of this "zeroth approximation" as the basic estimate, because in most of moduli 
space the curvatures are proportional to the metric up to 0(1) coefficients. In any case, 
this is the simplest distribution one can suggest which uses no data other than the effective 
theory itself, and is thus the "null hypothesis" in this class of problem. 

In comparing our actual results to this, perhaps the most striking difference is the 
growth of vacuum density in regions of large curvature, for example near conifold points. 
The simple physical explanation in this case is that we know that the structure of the 
potential, or equivalently the dual gauge theory description of the physics, can produce a 
hierarchically small scale A. Since the average spacing between vacua is expected to be A, 
we can fit more vacua into such a region. 

To make this quantitative, we need to understand the physics (or math) which makes 
the number of vacua finite. Here it is the tadpole condition, schematically AB = L where 
A, B are two conjugate types of flux, in dual terms controlling (say) the rank of the gauge 
group and the gauge coupling. A constraint AB = L translates a distribution dAdB into 
a scale invariant distribution L dA/A. 

In the large complex structure limit, the superpotential goes schematically as A + Bt^ 
for some power k = 1,2,3. This leads to vacua at A/B, and the distribution 

dA/A translates into d^r/Imr^, taking into account supersymmetry which brings in the 
complexified modulus. Such a measure is naively scale invariant, and would lead to a 
diverging number of vacua. But, it is important that the integration region for Rer is 
[0, 1), so in fact the distribution falls off for large Imr and is integrable; the total number 
of vacua with |r| > T goes as 1/T. One can also understand this as a consequence of 
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duality acting on the other fluxes, which leaves finitely many inequivalent choices (as in 
the toy example worked out in |j|]). 

Near a conifold, while the flux distribution is again roughly L dA/A, the dual gauge 
theory-type superpotential W = Az + Bzlogz leads to vacua at logz ~ A/B which have 
a jdloglogzp ~ d?z/\zlogz\'^ distribution. In fact, the distribution is rather similar to the 
previous one, with the identiflcation r = — logz. Whether this has deep signiflcance is not 
clear to us; the large complex structure limit and conifold limit are not dual and in other 
ways are not similar. But it means that the number of vacua N with a small scale z < e~'^, 
has a similar slow falloff, ~ 1 /T. 

In some sense, many vacua are "Kahler stabilized" - if the Kahler potential were 
different, they would not exist. An illustration can be found in flgure 4. As usual, this is 
less true for small e'^lW]'^. 

In section 5, we discussed the nature of the finite L distribution of vacua. Because 
of flux quantization, this is far more complicated than the smooth distributions which we 
have been computing. On the other hand, one can say a lot about it using the smooth 
distributions and the methods we discussed. The basic idea is to think of the smooth 
distribution as computing a the volume of a region in flux space which supports vacua. 
Physical flux vacua satisfying flux quantization are lattice points in this region. While 
the volume of the region controls the large L asymptotics, its other characteristics and in 
particular its surface area control the corrections to these asymptotics. 

In particular, the smooth distribution will always approximate the flnite L distribution 
for sufficiently large L. However, the mimimal L for this to be true depends on the size 
of the region in moduli space in which one is counting vacua, and can be estimated as 
L ~ K/r'^ for a ball of radius r. Thus, the number of vacua in a small region can deviate 
from the large L distribution and even the L^/^ asymptotic for the total number of vacua, 
until L becomes quite large. 

Another way to say this is that at finite L, to get a smooth density, one must average 
the number of vacua over regions in moduli space of size r ~ K/L. The actual number 
of vacua in smaller regions will show oscillations, with most regions having many fewer 
vacua, and a few having many more vacua to make up for this. We discussed the simplest 
example in section 5, and a one complex modulus example is discussed in |24]. 

It is not inconceivable that these considerations could significantly decrease the total 
number of vacua in interesting examples with many cycles, if the geometric factors such 
as total volumes of Calabi-Yau moduli spaces were sufficiently small. Work on computing 
these volumes is in progress. 

Let us move on to consider the results for specific types of vacua. First, we give the "null 
hypotheses" or simplest pictures which one might expect. First, F-type nonsupersymmetric 
vacua should be comparable in number to the super symmetric vacua (given our definitions, 
this is obvious for the D-type), up to factors like 2" where n is the number of moduli. A 
heuristic argument for this is that the potential V{z) is quadratic in W{z), so if W had 
been a polynomial of degree d, leading to {d— 1)" vacua, then V would be a polynomial of 
degree 2d, suggesting (2d — 1)" vacua. On refiection, the main problem with this argument 
is that the equations V = are real equations which typically have fewer than {2d — 1)" 
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solutions. Work on zeroes of randomly distributed polynomials with natural geometric 
distributions suggests that in some cases, most zeroes are real, while in others there 
are many fewer real zeroes, so this is inconclusive. Finally, the metastability condition 
might be expected to be weak, in the sense that if the masses of bosons are distributed 
symmetrically about zero, the the fraction of tachyon free vacua would be 2~^"'. 

In general, these expectations do hold, and the number of non-supersymmetric vacua 
is comparable to the supersymmetric vacua. However, the number of metastable nonsuper- 
symmetric vacua falls off drastically near conifold points, or more generally if the curvature 
on moduli space becomes large. This is rather surprising as naively the curvature contri- 
bution to the mass matrix Eq. (|1.8| ) goes the other way; the conifold point has positive 
curvature which raises the masses. The explanation for the D type vacua is not compli- 
cated; it has to do with the special form of the mass matrix which leads to a "mixing" with 
the dilaton-axion which forces a mode tachyonic, as explained in 3.1. Some experimenting 
with multimodulus models suggests that this phenomenon is specific to one modulus and 
the detailed form of this potential, again in the D breaking case. We observed the same 
phenomenon for the F breaking; it may have a similar explanation. 

We found in section 3.3 that the distribution of values of e'^jW^p or "AdS cosmological 
constants" is typically uniform near zero, so that the fraction of vacua with |A| < A* 
behaves as A^/LTs. A simple argument for this was mentioned in |16]: for a given value 
of moduli, the magnitude of is a single direction in "flux space," so the condition that 
it be small can be accomplished by one tuning of fluxes. 

However, the fact that this tuning can be made independent of the choice of moduli 
is non-trivial. In problems such as the attractor mechanism with only one type of flux or 
charge for each cycle, it is not true. The heterotic string also has one type of flux per cycle 
and is formally very similar; one cannot get small AdS cosmological constant just from 
fluxes in this case. Of course, there are many more variables having to do with the gauge 
fields, which might make small cosmological constant possible in this problem. 

The simple argument also misses a good deal of structure in the distribution; a good 
example is the behavior near the conifold point displayed in figure 4. 

The heuristic argument that "nonsupersymmetric vacua are as common as supersym- 
metric" actually has a precise realization here, in the existence of vacua with W = D^W = 
0. Since W = these are probably of limited physical interest, in particular it is hard 
to see how to stabilize their volume moduli. Other nonsupersymmetric vacua exist and 
are also roughly comparable in number. Strangely enough, the conditions of metastability 
and positive effective cosmological constant were incompatible, at least in the conifold and 
large structure limits. If this were generally true, and true in multi-modulus models, it 
would appear quite important. 

Of course, all this was only a first cut at an actual counting of vacua. At this point 
it seems quite possible that many of the nonsupersymmetric vacua which are stable under 
variations of the moduli, have more subtle instabilities or inconsistencies. This would 
obviously be important to understand. 

The scenario in which the numbers we are computing would be most meaningful would 
be one in which such effects drastically reduced the number of vacua, but in a way which 
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was more or less uncorrelated to the distributions we studied. If we grant this, certain 
physical expectations would result. First, one finds that a high scale of supersymmetry 
breaking is favored (of course, this could be the scale in a hidden sector). This is already 
true in the D breaking vacua, for which at best the scale Mf^^y is uniformly distributed, and 
in some cases [e.g. near conifold points) even disfavors small values. While we did not find 
F breaking vacua with positive cosmological constant, this may be an artifact of the one 
modulus case (or of the special limits we considered) . But a high susy breaking scale would 
be even more favored for F breaking vacua in models with many moduli, as one expects 
the many supersymmetry breaking parameters Fi = DiW to be roughly independent and 
uniformly distributed. 

The expectation that a high scale of supersymmetry breaking would lead to problems 
in tuning away the cosmological constant, favoring a low scale of breaking, is probably 
not true for these vacua. The point is that the contribution to the energy which does 
this tuning, the — Se'^^jVFp term or "AdS cosmological constant," gets contributions from 
many sectors of the theory, including sectors with unbroken supersymmetry, and has no a 
priori correlation to the supersymmetry breaking scale. Given a particular supersymmetry 
breaking scale, this contribution is then set by the requirement that A ~ 0, but since in 
general vacua exist with a uniform distribution of this parameter, this condition in itself 
does not favor any particular scale. 

Thus, these distributions seem to favor "moduli dominated gravity mediated super- 
symmetry breaking," in the terminology of One has a general prediction of a large 
gravitino mass, and risks the common problems of supergravity mediated breaking, such 
as non-universal soft breaking terms. Perhaps this is less of a problem in the D breaking 
models. 

Finally, the moduli space measure (the metric on the space of metrics) heavily disfavors 
large complex structure and large volume in the many modulus models, as we will discuss 
in more depth in ||TT| . 
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